Regular-to-Chaotic Tunneling Rates: Prom the Quantum to the Semiclassical Regime 



O 

(N 



U 



> 
m 

OO 

o 

ON 

o 



StefTcn Lock,^ Arnd Backer.^ Roland Kctzmcrick.^ and Peter Schlagheck^' ^ 

^ Institut fiir Theoretische Physik, Technische Universitdt Dresden, 01062 Dresden, Germany 
^Institut fiir Theoretische Physik, Universitdt Regensburg, 93053 Regensburg, Germany 
'Division of Mathematical Physics, Lund Institute of Technology, PBox 118, 22100 Lund, Sweden 

(Dated: March 18, 2010) 

We derive a prediction of dynamical tunneling rates from regular to chaotic phase-space regions 
combining the direct regular-to-chaotic tunneling mechanism in the quantum regime with an im- 
proved resonance-assisted tunneling theory in the semiclassical regime. We give a qualitative recipe 
for identifying the relevance of nonlinear resonances in a given ft-regime. For systems with one or 
multiple dominant resonances we find excellent agreement to numerics. 
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In mixed regular-chaotic systems the quantitative un- 
derstanding of dynamical tunneling, which refers to clas- 
sically forbidden transitions between phase-space regions 
that are separated by dynamical barriers [l|, represents 
one of the most challenging open problems in semiclas- 
sical physics. In the early 1990s it was found that 
tunneling rates between phase-space regions of regular 
motion are substantially enhanced by the presence of 
chaotic motion, which was termed chaos-assisted tunnel- 
ing. Such dynamical tunneling processes are ubiquitous 
in molecular physics [l| and were realized with cold atoms 
in periodically modulated optical lattices Q . 

The fundamental process in this context is tunneling 
between states in a regular island and the surrounding 
chaotic region. In the quantum regime, h < where 
Planck's constant h is smaller but comparable to the area 
of the regular island A, this process is dominated by a 
direct regular-to-chaotic tunneling mechanism 04^- An 
approach to determine these direct tunneling rates was 
recently given in Ref. [l^. It relies on a fictitious inte- 
grable system that resembles the regular dynamics within 
the island under consideration, and has been applied to 
quantum maps IC | , billiard systems [llj , and the annu- 
lar microcavity |12| . In the semiclassical regime, h <^ A^ 
however, the direct tunneling contribution alone is inca- 
pable to describe the observed tunneling rates. 

The existence of nonlinear resonances inside the reg- 
ular island leads to peaks and plateaus in the tunneling 
rates. These can be qualitatively predicted by the the- 
ory of resonance-assisted tunneling p^ . [l^ as shown for 
quantum maps [l^ , periodically driven systems [isl . [l6j , 
quantum accelerator modes , and for multidimensional 
molecular systems [l3j[3- Quantitatively, however, de- 
viations of several orders of magnitude appear (see Fig. [TJ 
gray line). Hence, both mechanisms alone give no reli- 
able quantitative prediction of tunneling rates in generic 
systems. 

In this Letter we derive a unified framework that com- 
bines the direct regular-to-chaotic and the resonance- 
assisted tunneling mechanisms. Beyond the direct tun- 
neling contribution there are additional contributions: 
They consist of resonance-assisted tunneling steps within 
the regular island and, in contrast to previous studies 



[1, 13 ' ^ final direct tunneling step to the chaotic sea; 
see insets in Fig. [TJ This leads to an excellent prediction 
of tunneling rates from the quantum to the semiclassi- 
cal regime. In particular, it includes the ft-regime, which 
is relevant for the search of experimental signatures of 
resonance-assisted tunneling. 

Specifically, we consider one-dimensional kicked sys- 
tems, described by the classical Hamiltonian H{p^ q, t) — 
T{p) + V{q) J2kez ^i^ ~ ^^'^ by the quantum time- 
evolution operator U = exp[— iV(g)/fi] exp[— ir(p)/S], 
which exhibit one major regular island embedded in the 
chaotic sea. This mixed phase-space structure gives rise 
to eigenstates of U that are mainly regular or chaotic, 
i.e. concentrated on a torus inside the regular region or 
spread out over the chaotic sea. However, they do have 
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FIG. 1: (Color online) Dynamical tunneling rates 7 from the 
innermost quantized torus (m = 0) of a regular island with 
one dominant 5:1 resonance. Numerical data (dots) is com- 
pared to the prediction of Eq. ((T)) (solid line) and previous 
results in the quantum regime due to the direct tunneling 
process (dashed line, Eq. @, Ref. [l3|) and in the semiclassi- 
cal regime due to resonance-assisted tunneling (gray solid line, 
Ref. UM ). The lower inset shows the phase space of the sys- 
tem [131 . The upper insets illustrate the dominant tunneling 
steps for three values of h. 
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small components in the other region of phase space. 
Consequently, a wave packet started on the mth quan- 
tized torus (m — 0,1,..., mmax ^ 1) of the island will, 
in the presence of absorbing boundary conditions in the 



chaotic sea, leak out from the island 



with a 



characteristic tunneling rate 7™. It describes the mean 
coupling to the chaotic sea, where fluctuations from in- 
dividual chaotic states are averaged. 

We now consider the presence of a prominent nonlin- 
ear r:s resonance (s oscillations match r driving periods) 
within the island. For small h, this resonance gives rise 
to couplings between different regular states [13, which 
compete with direct tunneling and lead to additional 
pathways into the chaotic sea. We derive (see below) 
the tunneling rate of the mth quantized state as 

7m — 7m ^ l^m.r I Im+r ^ \^m,r ^m,2r\ lm+2r ^ ■ ■ ■ ■ 

(1) 

The first term describes the direct tunneling process from 
the mth quantized torus to the chaotic sea with the rate 
7^, neglecting any influence from nonlinear resonances. 



In the second term the coefficient A> 



{r:s) 



describes one 



resonance-assisted tunneling step from the mth to the 
(m-|-r)th quantized torus via an r:s resonance, while the 
factor 7^^.^ describes the subsequent direct tunneling 
into the chaotic sea. The last term accounts for two 
resonance-assisted and one direct tunneling step. The 
three terms are visualized as insets in Fig. [T] in the h- 
regime where they are most relevant. Note, that the fi- 
nal step to the chaotic sea occurs due to direct tunneling, 
which is in contrast to previous studies 0, , where an- 
other resonance-assisted step is used. We now discuss (i) 

the direct tunneling rates 7JJ, 
and (iii) the derivation of Eq. ([IJ . 

(i) The direct tunneling rates can be predicted by an 
approach using a fictitious integrable system which 
has to be chosen such that its classical dynamics re- 
sembles the dynamics of the mixed system within the 
regular island as closely as possible. The eigenstates 
|m)o of the fictitious integrable Hamiltonian arc lo- 
calized in the regular region of H and decay outside. Us- 
ing Fermi's golden rule, 7^ = X^ch k|m)o,chP, in which 
= (V'ch|C^|m)o couples |m)o to different chaotic 
states, the direct tunneling rates arc given as [l^l 



(ii) the coefficients Am'P ^ 



7,t-||i'ch(C/-C/rcg)|m)o|| 



(2) 



where C/rog = exp(— ii7o/ft) and Pch is a projector onto 
the chaotic region. Here we assume that there are no 
additional phase-space structures within the chaotic sea 
that affect the direct tunneling rates. Note, that this 
approach is applicable for general systems with a mixed 
phase space [ll|, [l3|. Different methods for the deter- 
mination of the fictitious integrable system can be em- 
fused in 



ployed, based on the Lie transformation [21 



(used 



Figs.djanddl) or on the frequency map analysis 122 
in Fig. Details will be presented elsewhere [20|. 

(ii) The coefficients ^^^'^^^ in Eq. ([T|) depend on non- 
linear resonances. The classical dynamics of such an r:s 



resonance is described in the corotating frame by the ef- 
fective pendulum Hamiltonian [ij, [HI 



I,..s) + 2VriI)cos{re) (3) 



with rir:s = 2TTs/r, whcrc 2Vr(/) cos (r6') is the per- 
turbation in terms of the local action-angle variables 
(1,0). The simplest approach is the direct quantization 
of Eq. ([3|) in action-angle space by neglecting the ac- 
tion dependence of the coupling using Vns = V, (/r:s) 
and by making a quadratic approximation of Hq{I) 
around the action /].:s of the r:s resonance, leading to 
i/o(/) — fii :s(/ — /r:s) ~ (I ~ Ir.s)^ /^nir-.s- It is then found 
that the r:s resonance couples the mth excited state |m)o 
of the resonance- free island to the state \m + r)o with the 
coefficient 



A. 



{r:s) 
m^kr 



E„ 



(4) 



where E„i = [Im - /r:s]^/(2mr:s) + ^r-.slm with /,„ = 
h{m + 1/2) are the eigenvalues of the approximation of 
Hq. /r:s, Vr:s, and m,:;s Can be extracted from the classical 
phase space [13, [H] . 

(iii) To derive Eq. ([1} we start from Eq. ^ and apply 
quantum perturbation theory to obtain the perturbed 
regular states 



\m)o - 



-r)o - 



AtfA^Zl\m + 2r)o- 



(5) 

These are now inserted into Fermi's golden rule 7™ = 
J2ch b|m>,chP with v\m),ch = {Ah\U\m). Using the defi- 
nition of 7jjj and neglecting "off-diagonal" contributions 
of mixed terms (which at most amount to corrections of 
the order of a factor two) leads to Eq. ([T]). 

In the following we discuss two important improve- 
ments to the coefficients A^^^'^J^ in (ii): (iia) On one hand, 
the quadratic approximation of Ho{I) around / = Ir.s 
fails to predict the correct unperturbed energies Ejji es- 
pecially for resonances close to the border of the regu- 
lar island such that the peaks in the tunneling rates are 
located at a wrong position (see Figs. [5] and |31 dotted 
lines). Here, we use a semiclassical procedure to deter- 
mine a global approximation of the unperturbed system 
Ho{I) [iOl- Using its eigenenergies Em in Eq. Q leads 
to improved peak positions in the tunneling rates. 

(lib) On the other hand, the action dependence of the 
perturbation Vr{I) needs to be properly accoimtcd for. 
Here we use the fact that the Hamiltonian H{p,q,t) is 
analytic in the "harmonic oscillator" variables (P, Q) = 
\/2J sin 9 , cos 6) which result from {p,q) through an 
analytical canonical transformation. This yields Vr{I) oc 
1^ in lowest order in / and V,.(/)e^*''^ oc (Q ± iPy . 
Quantization in Q-space then amounts to replacing (Q + 
\P)/\/2h by the ladder operator a satisfying a|m)o = 
/m|m — 1)0. We therefore obtain in leading order 



(6) 
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with m' = m + {k— l)r. These new couphng parameters 
V™g''^^ depend on H and on the quantum number m of the 
coupled regular state. They replace Vi:s in Eq. ^ and 
improve the tunneling rates in between the peaks. 

We now predict dynamical tunneling rates with Eq. ([T]) 
using (i) the direct tunneling rates -f'^ from Eq. ([2]) de- 
scribing the decay of the regular states by ignoring the 
influence of resonances, (iia) the improved eigenenergies 
Em of Ho{I) and (lib) the action dependent couplings 
■pm.fcr fj-Qjjj gq 1^ Fig. [T] shows a comparison between 
this theoretical prediction (solid line) and numerical data 
(dots), obtained using absorbing boundary conditions in 
the chaotic sea. In contrast to the sole application of 
direct tunneling with Eq. ([2|) (dashed line) or resonance- 
assisted tunneling according to Ref. (gray line) we 
find excellent agreement from the quantum to the semi- 
classical regime. The phase space of the system is de- 
picted in the inset of Fig. [TJ It shows one dominant 5:1 
resonance inside the regular island. All other resonances 
arc small and do not contribute to the tunneling rates for 
the shown values of h. We find that the direct tunnel- 
ing rates following from Eq. ^ explain the numerical 
data for 1/h < 20. For larger 1/h the 5:1 resonance is 
important and gives rise to two characteristic peaks cor- 
responding to the coupling of the ground state m = 
to m = 5 and m = 10. Our theoretical prediction excel- 
lently reproduces the peak positions and heights. 

As a next example we consider a situation where three 
dominant resonances exist, namely a 5:1 resonance near 
the center of the regular island, an 11:2, and an outer 
6:1 resonance, see the inset in Fig. [2l Here, multireso- 
nance couplings occur and the tunneling rates from the 
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FIG. 2: (Color online) Dynamical tunneling rates 7 from the 
quantized torus m = 2 of a regular island with three domi- 
nant resonances. Numerical data (dots) is compared to the 
prediction of Eq. ([1} (solid line). For comparison the direct 
tunneling contribution (dashed line, Eq. ([2]), Ref. [l3]) and 
the result of Eq. ([T]) without the improvements (iia) and (iib) 
is presented (dotted line). The inset shows the regular is- 
land of the system [l^ and the dominant 5:1, 11:2, and 6:1 
resonances which cause tunneling steps as indicated by labels. 



torus TO = 2 are determined with Eq. ([l} by a summation 
over all relevant r:s resonances. For example the peak at 
l//i ss 84 in the tunneling rates is caused by the coupling 
of the state to = 2 to the 18th excited state via the 5:1 
and the 11:2 resonance. In Fig. [2] we compare numeri- 
cal data (dots) to the prediction of Eq. ([T]) and find very 
good agreement. For comparison we also show the rates 
(dotted line) that result from using Eq. ([T]) together with 
Eq. ([4]) without the improvements (iia) and (iib). We 
still find overall agreement to the average decrease of the 
tunneling rates with 1/h. However, the reproduction of 
the positions and heights of the individual peaks requires 
to use a global approximation for the unperturbed ener- 
gies (iia) and to take into account the action dependence 
of the coupling matrix elements (iib). 

The paradigmatic model of quantum chaos is the stan- 
dard map [T{p) = pV2, V{q) = -K/{iTT^) cos{2Trq)]. 
For K = 3.5 it has a large generic regular island with a 
dominant 6:2 resonance and further relevant 14:5 and 8:3 
resonances. These are the largest low-order resonances. 
Fig. [3] shows the comparison between numerical data and 
the prediction of Eq. ([1]) far into the semiclassical limit. 
We find quantitative agreement over 30 orders of mag- 
nitude in 7. The direct tunneling process is relevant for 
1/^ ^ 30 beyond which several regimes for the tunnel- 
ing rates caused by different resonances arc identified. 
We attribute the small deviations below the prediction 
of Eq. Il]) to destructive interference of different tunnel- 
ing sequences leading to the same final state. 

An important question is when the resonances become 
relevant, i.e. where is the transition between the direct 
tunneling regime h < A, where resonances play no role. 
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FIG. 3: (Color online) Dynamical tunneling rates 7 from the 
innermost quantized torus m = of the regular island for the 
standard map at K = 3.5. Numerical data (dots) is compared 
to the prediction of Eq. ([T]) (solid line). For comparison the 
direct tunneling contribution (dashed line, Eq. ((2|, Ref. [TO]) 
and the result of Eq. ^ without the improvements (iia) and 
(iib) is presented (dotted line). The inset shows the island 
with the dominant 6:2, 14:5, and 8:3 resonances and the ar- 
rows indicate the regimes where they start to become relevant. 
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and the resonance-assisted tunneling regime /i <C -A. The 
latter requires that the (m+r)tli quantized torus exists 
inside the island, i.e. 



A 

— > m 
h 



1 



(7) 



The position /ipoak of the first peak in jm arises when 



Eq. (HI diverges for ~ E., 



rhQy-s = 0. Using the 



quadratic approximation for Ho{I), we find 



^peak 



(8) 



with ^1:8 = 27r/r:s- Thc first peak appears when the r:s 
resonance encloses m+r/2 quantized tori. However, the 
influence of the resonance may appear much earlier than 
f//ipeak. We define the transition point l//iros when the 
first two terms in Eq. ([I} are equal. Using V2TOr:sVr:s = 
A^i.:s/16 where A^i:s is the area enclosed between 
the separatrices of thc r:s resonance, we obtain 



1 



rhr 



•Ar:s y '7m(^rcs) (^rcs/^pcak 1) 



128 

7^2 



(9) 

This criterion explicitly involves the size AAr-.s of the 
r:s resonance chain measured with respect to r Planck 
cells and to the area Ar-.s enclosed by the resonance, 
as well as the ratio of the direct tunneling rates from 
the (m+r)th and the mth quantized torus. As thc left- 
hand side of Eq. ([9]) is expected to display a monotonous 



increase with l//ircs towards the singularity at l//ipeak, 
resonance chains with a large area A^i.:s will lead to a 
lower transition point 1 / h^cs at which the crossover from 
direct to resonance-assisted tunneling appears. While, 
e.g., the first peak of the 5:1 resonance in Fig. [2] appears 
at l//ipcak ~ 130 it dominates the tunneling process al- 
ready at l//ircs ~ 40 which is even before the contribu- 
tions from other resonances set in. Equations ([5]) and 
([9]) confirm the intuition that low-order resonances with 
small r,s and large A^i-:s are most relevant. 



In conclusion, we have combined an improved 
resonance-assisted tunneling theory with the theory of di- 
rect tunneling using a fictitious intcgrable system to pre- 
dict dynamical tunneling rates from the quantum regime, 

^ A, to the semiclassical regime, h <^ A. Excellent 
quantitative agreement with numerically determined tun- 
neling rates is found on the level of individual resonance 
peaks, which emphasizes the validity of the underlying 
direct and resonance-assisted mechanisms. We therefore 
expect that these mechanisms leave their characteristic 
traces in semiclassical approaches based on complex clas- 
sical trajectories [l^l, and allow one to understand and 
predict tunneling rates also in more complex, multidi- 
mensional quantum systems. 
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